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ABSTRACT 

IRAS04325+2402 (heraftcr IRAS04325) is a complex protostellar system hosting two 
young stellar objects (AB and C in the following) at a separation of 1250 AU. Here we 
present new deep Gemini imaging and spectroscopy for the system covering the wave- 
length regime from 1-12 /im as well as Sub-Millimeter Array interferometry at 870 /xm, 
in combination with Spitzer and literature data. Based on this rich dataset we provide 
a comprehensive picture of IRAS04325 over scales from a few AU to several parsec. 
Object AB is a low-mass star with a disk/envelope system and an outflow cavity, which 
is prominently seen in infrared images. Object C, previously suspected to be a brown 
dwarf, is likely a very low mass star, with an effective temperature of ^ 3400 K. It 
features an edge-on disk and an elongated envelope, and shows strong indications for 
accretion and ejection activity. Both objects are likely to drive parsec-scale molecular 
outflows. The two objects are embedded in an isolated, dense molecular cloud core. 
High extinction, lack of X-ray emission, and relatively high bolometric luminosity ar- 
gue for a very young age below 1 Myr. The disk/outflow systems of AB and C are 
misaligned by ^ 60deg against each other and by 80 and 40deg against the orbital 
plane of the binary. The system might be a good case for primordial misalignment, 
as opposed to misalignment caused by dynamical interactions, because the outflow 
direction is constant and the realignment timescale is likely larger than the system 
age. This favours turbulent fragmentation, rather than rotational fragmentation, as 
the formation scenario. We show that the spectral energy distributions and images for 
the two objects can be reproduced with radiative transfer models for disk/envelope 
systems. Our analysis provides reassurance in the established paradigm for the struc- 
ture and early evolution of YSOs, but stresses the importance of developing 3D models 
with sophisticated dust chemistry. 

Key words: stars: low-mass, brown dwarfs, stars: circumstellar matter, stars: pre- 
main sequence, stars: formation 



1 INTRODUCTION 

The study of protostars and their environment requires the 
analysis of radiation across a broad range of wavelengths, 
including the infrared and the submm/mm regime. At the 
same time, it is desirable to achieve high spatial resolution 
to characteriz e the small-scale stru ctures close to the sources 
(see review by I Watson et al.ll2007l ). However, only a few ob- 
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jects have been studied at uniform and high resolution over 
the full wavelength domain of relevance (1-1000 fj,m). Testing 
the current paradigms for collapse, accretion/outflow, and 
disk evolution with high-resolution multi-wavelength stud- 
ies of individual objects is thus of major importance. 

In this paper we present a case study for the protostellar 
system IRAS04325-I-2402 (hereafter: IRAS04325), shown in 
Fig-Ill Origina lly identified as LDN 1535 in the catalogue of 
dark clouds bv lLvndsl (|l962r ). it has been studied extensively 
over three decades (117 publications according to SIMBAD). 
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Figure 1. K-band close-up view of IRAS04325 obtained with 
Gemini-North and NIRI. The image is shown in two different 
intensity scales: the left panel shows the point sources and bright 
extended emission close to the southern source AB, and the right 
panel illustrates the faint structure between AB and the northern 
source C. 



A number of publications starting in the 1980s has focused 
on the molec ular gas in L1535 and the large-scale envi- 
ronm ent (e.g. lUngerechts et al.ll 19821 : iGoldsmith fc SernvakI 
[l98J). After the identification of an IRAS source associ- 
ated with this cloud (Parker 1988, 1991|), the object has 
been analysed in the infrared domain (e .g. Hodapd Il994l'l 
and in the submm/mm continuum (e.g . iLadd et al.l Il99ll : 
lOhashi et al.l Il996l : ISaraceno et al.l 119961 ) , tracing the dust 
in the system. IRAS04325 was found to harbour a protostar 
with subsolar lumino sity and a disk seen at high inclination 
IjKenvon et al.lll993al lbr) , embedded in a cloud core with a di- 
ameter of several tens of thousands of AU. In ad dition, large- 
scale outflow activity has been identified (e.g. iHever et al.l 
ll987l : lMoriartv-Schieven et al.l ll992': Wang et al. 2001). The 
source is established as a prototypical small-scale site of star 
formation. 

Until 1999, the structures within the IRAS beam have 
not been resolved, which results in various problems. The 
near-infrared images clearly show a bright nebulosity, in- 
dicating that the fluxes from IRAS and single-dish obser- 
vations are contaminated by extended emission. As a result 
the properties of the central source and its close environment 
are poorly constrained. Confusingly, the direction of one of 
the large-scale outflows does not agree with the orienta- 
tion of the infrared nebulosity, int e rpreted as outflow cavity, 
as po inted out by iKenvon et al.l (|l993al ). iHartmann et al.l 
([1999|) presented HST imaging of the object, revealing com- 
plex substructure. In particular, the images show the pres- 
ence of a second well-separated point source (IRAS04325C), 
which has a resolved disk seen edge-on. The main light 
source IRAS04325AB is speculated to be a binary. In ad- 
dition there are various scattered light features. 

Originally fueled by the suspicion the second source 
IRAS04325C might be a proto brown dwarf, we have con- 
ducted our own obs ervations of the system. In our flrst paper 
l|Scholz et al.ll2008h . we detect the second component with 
1.3 mm interferometry and show the first infrared spectrum 



for this object. Here we present a comprehensive dataset for 
the full IRAS04325 system, including new near- and mid- 
infrared imaging, near-infrared spectroscopy for both com- 
ponents, and submillimeter interferometry. We consistently 
assu me a distance of 150 p c for the Taurus star forming re- 
gion (JLoinard et al.l 120071 ) when converting from apparent 
to physical separations. All physical distances are given in 
projection against the plane of the sky and should be seen 
as lower limits to the actual distances. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Near-infrared imaging 

Near-infrared images in J-, H-, and K'-band were taken on 
the 2 3 August 2008 with the NIRI camera l|Hodapp et all 
|2003| ) at the Gemini-North 8.1 m telescopfE]. All Gemini ob- 
servations for IRAS04325 were taken in queue mode in the 
band-2 program GN-2008B-Q-73. We used the f/6 camera 
with a scale of O'.'12/pix and 120" x 120" FOV. A 7-step 
dither pattern was carried out in each band with 15" offsets 
and 10 sec integration time per position. The total on-source 
time in each band was 70 sec. The seeing measured from 1-3 
well-detected point sources outside IRAS04325 was 0'.'5. 

The reduction of the NIRI J-, H- and K'-band images 
proceeded as follows: a sky image for each filter was created 
using the dithered science images. The sky was then sub- 
tracted to each science image using the ESO eclipse package 
l|Devillardlll997l ). We carried out a flatfield correction using 
domeflats. The alignment and coaddition of th e images was 
done with the pro grams Scamp (jBertirJ 120061 ) and Swarp 
IjBertin et al.ll2002h to produce the final science stacks. 

In the same Gemini program, we obtained L'- (3.8 ^m) 
and M'-band (4.7 /xm) images with NIRI's f/32 camera 
(O'.'02/pix) and a FOV of 22.4" x 22.4". In the L-band a five 
step dither pattern with 3" offsets was executed twice with 
20 coadded 1 sec exposures at each position, which gives a 
total on-source time of 200 sec. The same pattern was carried 
out 13 times in the M-band with 30 coadded 0.5 sec expo- 
sures, in total 25 min. These observations were obtained on 
the 3rd of December 2008, including a set of images for the 
bright AO star HD36719 for calibration purposes. All data 
was taken at low airmass (< 1.2); the airmass difference 
between science and calibration target was < 0.1. 

For the reduction of the L- and M-band data we followed 
the standard recipes based on the tools in the Gemini/NIRI 
package within IRAF. Flats and sky images were created 
based on the dithered science images, discarding the flrst 
frame. In the L-band all nine available images were used 
for the sky frame. The M-band sequence is long enough for 
the sky to change significantly. For each individual frame 
we used the nine frames taken closest in time to create 



^ The Gemini Observatory is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partnership: 
the National Science Foundation (United States), the Science and 
Technology Facilities Council (United Kingdom), the National 
Research Council (Canada), CONICYT (Chile), the Australian 
Research Council (Australia), Ministrio da Cincia e Tecnologia 
(Brazil) and Ministerio de Ciencia, Tecnologa e Innovacin Pro- 
ductiva (Argentina) 
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the sky image. After subtracting the sky the science frames 
were divided by the normalized flatfield. The resulting im- 
ages were mostly flat with the exception of the first and 
last two frames in the M-band suffering from imperfect sky 
subtraction. These images with uneven background were dis- 
carded. To coadd the individual images we used the standard 
combine routine in IRAF. The images of the standard star 
for L- and M-band were handled in exactly the same way. 



2.2 Mid-infrared imaging 

IRAS04325 was obs erved with Gemini-North and Michelle 
l|Brvson et al.l 11994 ) in the N'-band with a central wave- 
length of 11.2 /im on the 3rd of November 2008. To sample 
the sky properly the standard procedure including chopping 
and nodding was used. In short, the pointing was chopped 
at 4.2 Hz between the target and an adjacent position 15" 
to the east by moving the secondary mirror. In addition the 
whole telescope was nodded in the same direction in larger 
time intervals. The chop throw was chosen to avoid the ex- 
tended emission from IRAS04325 from the sky field. The 
exposure time per frame was 0.02 sec; in total we integrated 
2 X 590 sec on source. The airmass for the two science inte- 
grations was 1.07 and 1.02. In addition, the standard star 
HD20893 was observed at airmass 1.02. 

For the standard reduction we used the Gemini/MIDIR 
package in IRAF and followed the 'cookbook' for Michelle 
provided by Gemini. This includes the tasks miregister to 
prepare and reformat the data cubes and mistack to average 
by coadding the signal from each nod position and dividing 
by the number of frames. A residual pattern was removed 
by fitting each column with a constant and subtracting this 
fit. Source C is weakly visible only in one of the resulting 
stacked frames - the one taken at lower airmss. The peak 
intensity for object C is about three times the noise in the 
sky background. The signal does not improve significantly 
when coadding the two images. Only the image with detec- 
tion was used for the photometry. The standard star frames 
were reduced in the same fashion as the science images. 



2.3 Near-infrared spectroscopy 

Object AB was observed with SofI (|Morwoodlll997l ) in spec- 
troscopy mode at the ESO/NTT in La Silla. We obtained 
this spectrum and the associated calibration files from the 
ESO archive. The low-resolution Red Grism was used with 
coverage from 1.5 to 2.5 ^m and a 0'.'6 slit, which results 
in a resolution oi R ^ 1000. Two pairs of nodded integra- 
tions (ABBA) were observed with a nodding step of 20" 
and 225 sec exposure time per position. The total on-source 
time was 900 sec. A telluric standard star (HIP 23719, spec- 
tral type GIV) was observed immediately after the science 
target. 

We carried out a standard reduction, including subtrac- 
tion of the nodded pairs (A-B) to remove the sky and fiat- 
field correction. The four spectra were extracted with a 1'.'2 
aperture, wavelength calibrated using Xenon arc spectra, 
and coadded. The telluric spectrum was treated in the same 
way. The combined spectrum for IRAS04325AB was then 
divided by the telluric spectrum and multiplied by the solar 
spectrum. The final science spectrum is shown in Fig. [4] 



We obtained a spectrum for object C on Oct 14/15 
2008 with Gemini-North and NIRI, using the f/6 camera, 
the 6-pix slit mask (slit width 0'.'72) and the H- and K-band 
grisms. This gave a resolution of i? ~ 500 over a wavelength 
range of 1.4-2.5 /xm. The slit was aligned in E-W direction. 
The total on-source time of 3000 sec per band was split in 10 
single integrations. To correct for the variable background 
emission, the object was moved along the slit in steps of 
2" between the single exposures. For telluric correction we 
observed the early A stars HIP17791 and HIP23088 in the 
same observing nights. 

The NIRI spectra were extracted and calibrated using 
our own IDL procedures. The frames were flat flelded using 
an internal fiat taken immediately after the science frames. 
Consecutive exposures are subtracted to remove the sky 
background. The 10 individual spectra are then extracted 
and median-combined. The same operation is performed for 
the telluric calibration stars. From the standard spectra we 
derived a telluric spectrum by dividing through a black body 
spectrum with a temperature of 10000 K and interpolating 
over the hydrogen lines. The science spectrum was divided 
by the telluric spectrum. The wavelength calibration was 
carried out based on the bright OH lines which are visible 
in the (unreduced) science images. The final H- and K-band 
spectra are shown in Fig. [5] 

2.4 Archival images 

To complement our dataset, we used archival images for the 
Taurus region. Our target and its environment is covered 
by the TAP Survey, carried out with the near-infrared cam- 
era WFCAM at t he United Kingdom Infrared TelescopsQ 
JDavis et al.ll2008r ). The survey provides a database of near- 
infrared images in the K-band and the narrow-band fil- 
ter centered on the 1-0 S(l) hydrogen emission feature at 
2.122 /xm (in the following called: H2 band), which we use 
to look for shocked emission caused by molecular outflows. 
In addition, we downloaded the pipeline-reduced PBCD 
Spitzer images from the 3rd delivery (DR3) of the Tau- 
rus Spitzer Legacy project (PI: Deborah Padgett, program 
id #3584) from the NASA/IPAC Infrared Science Archive. 
This includes the IRAC data in the four channels at 3.6, 4.5, 
5.8 and 8.0 fim, as well as the MIPS images at 24 and 70 fim. 

2.5 Submillimeter interferometry 

A field centered on IRAS04325 was observed with the SMAQ 
on 20 October 2008 using eight antennas in a compact 
configuration that provided baseline lengths of 11 to 77 m. 
The weather was good, with system temperatures of 200 
to 250 K (DSB) near transit. The phase center was RA 
Q^hggmggsgg^ DEC +24°08'27'.'6 (J2000), which is the po- 
sition of IRAS04325C. The correlator was configured for 

^ The TAP survey is hosted by the Joint Astronomy Center. Im- 
age reduction was carried out by the Cambridge Astronomical 
Survey Unit for the WFCAM Science Archive. For more informa- 
tion, see 'http : //www . j ach . hawail . edu/UKIRT/TAP/ 
^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Figure 2. SMA 870 ^m dust continuum image showing the two 
sources AB and C and associated extended emission. The ellipse 
in the lower left corner shows the 2''0 X 1'.'6 PA — 68deg beam. 
The lowest contour and step are 4.4mjy/beam (twice the rms 
noise level). Negative contours are dashed. 



the full 2 GHz bandwidth, including one spectral chunk 
of 82 MHz in the middle of the upper sideband devoted 
to the CO J=3-2 line at 345.786 GHz with channel spac- 
ing 0.20MHz (0.176 km s~^). The SMA primary beam 
FHWM is 35" at this frequency. Calibration of complex 
gains was performed by interleaving 3 minute observations of 
the quasars 3C111 (4.30 Jy) and J0530+135 (0.84 Jy) with 
15 minute observations of IRAS04325 over the hour angle 
range —5 to +5. The flux scale was set using observations of 
Uranus, with an estimated uncertainty better than 20%. The 
passband was calibrated using observations of the strong 
sources 3C454.3 and Uranus. The data calibration was done 
with the MIR/IDL software, followed by standard imaging 
and deconvolution tasks within the Miriad package. For the 
continuum image, the beam size (robust=0) was 2'.'0 x l'.'6 
at PA — 68deg, and the rms noise 2.2mJy/beam, after com- 
bining the lower and upper sidebands (effective bandwidth 
~ 4 GHz). For the CO line images, the beam size (natu- 
ral weight) was 2'.'l x l'.'8 at PA — 56deg, with rms noise 
0.26 Jy/beam in each spectral channel. The continuum im- 
age is shown in Fig. (2] the CO line maps for a sequence of 
velocity bins are plotted in Fig. O Since the interferome- 
ter observations are sensitive only to small scale structure 
(< 10"), the sources seen in the CO maps are not isolated 
and should be seen as local overdensities and temperature 
enhancements. 



3 QUANTITATIVE ANALYSIS 

3.1 Object AB: NIR spectroscopy 

The spectrum of the source AB is red and entirely feature- 
less (see Fig. |4]), without the bro ad water absorption bands 
usually seen in late-type objects (jCushing et al.ll2005h . The 



slope from 1.5 to 2.3 ^m is well approximated by a straight 
line. There may be a hint of CO absorption bands at 
> 2.3 ^m. 

The lack of features precludes a detailed spectral anal- 
ysis. Instead we simply compare the slope with reddened 
spectra for late-type object s from the IRTF/SpeX Spectral 
Library (jRavner et al.ll2009l ). G-, K- and early M-type spec- 
tra reddenend with Ay ~ 30 provide a good match. As an 
example we show in Fig. U the spectrum for the K7V star 
HD237903. Towards mid-M spectral types the fit becomes 
progressively worse, no matter which extinction is chosen. 
In particular, the slope over H- and K-band is no longer 
constant, due to the water absorption feature between the 
two bands. In addition, the CO bands and water absorption 
cause the flux level to drop at > 2.2 /im. The mismatch is 
shown in Fig. [4] where we compare with the spectrum of 
the M4V star Gl 299 reddened with Av = 30 mag. Thus, a 
spectral type later than early M seems unlikely. 

This limit on the spectral type, however, needs confir- 
mation, particularly because the spectrum may be affected 
by additional continuum emission from the disk or accretion. 
The slit orientation in E-W direction avoids the scattered 
light features north and south of the source (see Fig. (Tjl , but 
veiling from the inner disk (within 0"5) cannot be excluded. 
Lacking high-resolution spectroscopy, we cannot put limits 
on the veiling. Typical excess continuu m in T T auri stars 
make s the near-infrared spectrum redder (IFolha fc EmersonI 
[l999|), i.e. the spectral type may be earlier than estimated 
or the extinction may be overestimated. 

3.2 Object C: NIR spectroscopy 

The spectrum of object C shows features of the photosphere, 
in particular a flattening at the red end of the H-band and 
a triangular shape in the K-band (Fig[5]). The spectra in H- 
and K-band were compared separately to model spectra. We 
neglect that the spectra might be affected by scattering and 
veiling. Both effects are expected to be smoothly and slowly 
varying functions of wavelength, thus they should not affect 
the spectrum in individual bands significantly. 

For the co mparison we used a series of AMES-DUSTY 
model spectra (jAUard et al.l 1200 Ih with effective tempera- 
tures ranging from 2000 to 4000 K and logg = 3.5 or 4.0, 
as expected for young stars. We varied Toff and Ay and 
aimed to obtain a consistent solution for the two bands. We 
find a decent match between observed and model spectra 
for Toff = 3400 ± 200 K and extinction Av = 11±2 mag. As 
seen in Fig. [S]this fits H and K bands equally well. 

It is impossible to match the observed spectrum with 
a model for an object close or bel ow the subst e llar b ound- 
ary (Tcft < 3000 K), as claimed bv lScholz et~ai] (|2008D . par- 
ticularly not in the H-band. Young brown dwarfs exhibit 
a characteristic tri angular shaped H-band spectrum due to 
wate r absorption IIKirkpatrick et al.l 120061 : iBrandeker et al.l 
120061 : IScholz et al.ll2009l ) which is not seen in our data. We 
conclude that IRAS04325C is unlikely to be a brown dwarf; 
its spectrum is more consistent with being a very low-mass 
star with an effe ctive temperature aro und 3400 K or spectral 
type of M3-M4 (JMentuch et al.ll2008l '). 

Various emission features associated with molecular hy- 
drogen are clearly detected in the K-band, particulary 1-0 
S(2) at 2.034, 1-0 S(l) at 2.122, 1-0 S(0) at 2.223, 1-0 Q(l) at 
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Figure 3. Contour maps of the continuum-subtracted CO line observations with SMA. The contours are in steps of 1.5 Jy/beam, negative 
contours are dashed. The LSR velocity in kms~^ is indicated for each panel. The positions of the compact sources AB and C, as measured 
in the continuum map, are marked with crosses. 



2.407, and 1-0 Q(3) at 2.424 /im. In addition, there is a Br 7 
emission feature at 2 .166 Atm, commonl y associated with on- 
going accretion (e.g. iNatta et al.ll2006l ). The H-band shows 
the [Fell] hnes at 1.644, 1.669, 1.677, and 1.748 /xm (the lat- 
ter might be blended with a hydrogen feature). While the 
lines of molecular hydrogen can originate in the disk or in 
a jet, the [Fell] emission can only occur in the jet. Thus, 
our spectra prove ongoing accretion and outflow activity in 
IRAS04325C. From the way we have analysed the spectra it 
is certain that the emission originates in an area within the 
PSF of the point source, i.e. within a radius of O'.'S (75 AU) 
around the source. 



The [Fell] lines p rovide a constraint on the electron den- 
sity W e in the jet (e.g. iNisini et al.ll2002l : iGarcia Lopez et al.l 
[2003)- The flux ratio between the well-detected lines at 
1.644 and 1.677 /xm is 2.9, which gives n^ > fO'^ cm"^. Note 
that the [Fell] line ratio saturates for larger densities. Con- 
sistently, the lower limits to the flux ratios of 1.544 vs. 
1.644 /im and 1.600 vs. 1.644 ^m yield densities > 0.5 • lO'^ 
and > 10*cm~^, respectively. Such high values are usually 
found at the base of the jet, while t he density typically 
drops further away from t he source (e.g. iTakami et al.ll2006l : 
iGarcia Lopez et al.l 120081 ). Thus the analysis confirms that 
the [Fell] emission originates in a jet close to IRAS04325C. 



The Br 7 line is used t o derive an e stima te for the 
mass accretion rate following iNatta et al.l (l200g). The ap- 
proximate equivalent width in this line is 7 ± 2 A. Scaling 
with the K-ba nd magnitude of 8. 5 mag for a 3400 K star 
at 1 Myr from iBar affc ct al.j (|l998l ) gives a line luminosity 
of logiLBi-y / Lq) ^ —4.1. With the empirical relations by 
ICalvet et al.l (|200J) this translates into an accretion lumi- 
nosit y of log(Lacc/I'©) ^ —0.8 (see also iMuzerolle et al.l 
1 19981 ). The mass accretion rate is then M = Li^ccR/GM. 
Assuming a mass of 0.4 M© and a radius of 1.5 R© results in 
an accretion rate in the range of (1.0 ± 0.3) • 10~* M© yr~^, 
comparable to accretion rates derived for T Tauri stars but 
highe r than those for young brown dwarfs l|Mohantv et al.l 

[iooi). 
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Figure 4. Spectrum for AB obtained with NTT/SOFI. Overplot- 
tcd with dashed lines are the spectra for a K7V star (above) and 
a M4V star (belo w the science spectrum), from the IRTF/SpeX 
Spectral Library IjRavner et al.l l2009.). Both comparison spectra 
are reddenened with Ay = 30 mag. 



3.3 Photometry of point sources 

We derived magnitudes/fluxes for the two components of 
IRAS04325 from the Gemini near- and mid-infrared images 
as well as from the submm continuum map. 

For the JHK images we adopted an aperture of l'.'2 
(10 pixel) and a background annulus between 1'.'2 and 1'.'8 
(10-15 pixel). The field contains one 2MASS source outside 
IRAS04325 with J = 17.98, H = 16.41, and K = 14.64 
which was used to shift the instrumental magnitudes into the 
2MASS system. The resulting values are listed in Table[T] In 
the J-band there is diffuse emission at the position of object 
C, but no clear point source; we consider the flux measure- 
ment to be an upper limit. The previously published pho- 
tomet ry of IRAS043 25 has been obtained with larger aper- 
tures IjHartmann et a l., 1999: Connolley ct al. 2008). Hence, 
these literature fluxes are significantly larger than ours be- 
cause they include extended emission. 

For the L'-, M'-, and N'-band photometry we used aper- 
tures of I'.'O or l'.'2 and a sky annulus similar to the one for 
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Figure 5. H- and K-band spe ctrum from Gemini/NIRI (solid line) in comparison writli a model (DUSTY, 3400 K, logg = 4.0, Ay = 
22 mag, see lAUard et al.l 1120011 )). The model (dashed line) is shifted in y-direction for clarity. The emission lines detected on the Za level 
are marked with dashdotted lines (see text for more information). 



Table 1. IR photometry for the point sources in IRAS04325. If 
not specified, the uncertainty in the fluxes is ~ 10%; this includes 
the calibration er ror. For more detail s see Sect. [331 The fluxes at 
1300 Atm are from lScholz et al.l ||2008| ). 



Band A (/^m) AB 



C 



unit aperture (") 



J 


1.25 


16.54 


< 18.6 


mag 


1.2 


H 


1.65 


14.37 


16.92 


mag 


1.2 


K' 


2.15 


11.87 


14.75 


mag 


1.2 


L' 


3.78 


10.36 


13.4 ±0.2 


mag 


1.0 


M' 


4.68 


9.37 


12.3 ±0.4 


mag 


1.0 


N' 


11.2 


65 


3.4 ±1.1 


mjy 


1.2 




870 


55 


41 


mjy 


3.0 




1300 


lOitf 


10±3 


mJy 


2.0 



JHK. In these bands the extended emission is less problem- 
atic, thus the choice of the sky annulus does not significantly 
affect the results. The L'/M'-band magnitudes were shifted 
to the standard system based on the values for the cali- 
bration star HD367 19, which has 5.4 m ag throughout this 
wavelength regime (|Leggett et al.M2003l ). For the N'-band 
magnitudes we used the star HD20893 for calibration, for 
which an N'-band flux of 4.3 Jy is reported on the Gemini 
websitqj, based on lCohen et al.l (|l999r ). 

The fluxes in the 870 /^m continuum map, measured 
within a 3'.'0 aperture, are 55 and 41 mJy, for AB and C 
respectively. 

The uncertainties were estimated by varying the size of 
the aperture and the sky annulus and checking the changes 
in the fluxes. Combined with the uncertainty in the fluxes 
used for the calibration, this mostly yields errors < 10%. 
The exceptions are the L'- and M'-band images for object 
C, where the errors are substantially larger. 

In Fig. [6] we analyse the NIR colours and magnitudes of 
the two sources (flUed symbols). In addition, we estimated 
the photospheric colours and magnitudes by comparing the 



^ |http : //www ■ gemlnl . edu/sclops/lnstruments/] 
mir/MIRStdFluxes . txt 
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Figure 6. Colour-magnitude diagram illustrating the interpreta- 
tion of the near-infrared photometry. Filled symbols show mea- 
sured photometry; empty symbols the approximate photospheric 
fluxes obtained by compari ng the constraints o n T^s with the 
1 Myr evolutionary track by iBaraffe et al.l 1 11998^ . Squares corre- 
spond to object AB and octagons to object C. Th e dotted line i s 
the reddening path for Av = 1-33 with the iMathij ill 9901 ) 
extinction law. Th e dashed line show t he standard scattering 
model MRN12 from lKenvon et al.l l ll993bl ) with inclinations rang- 
ing from cos i = 0.975 — 0.0. The solid line is the scattering model 
plus extinction. 



constraints on the e ffective temperatures with the 1 Myr evo- 
lutionary track by iBaraffe et al.l (|l998l ) (empty symbols). 
For source C, we used T^s = 3400 K, as derived by compar- 
ing the spectrum to models (Sect. I3.2|l . For object AB we 
only have a limit on the spect ral type (< early M) , which 
corresponds to T^s > 3600 K (JMentuch et al.ll2008l ). Thus, 
the plotted flux constitutes a lower limit as well. These pho- 
tospheric datapoints may be uncertain in the y-axis by up 
to ±1 mag due to the unknown age and the uncertainty in 
Tes , but this does not affect the following interpretation sig- 
nificantly. 

For object AB (square symbols) the colours are well ex- 
plained by the extinction path (dotted line) , indicating a line 
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of sight extinction of Ay ~ 30 mag, in agreement with our 
constraint from the spectrum (Sect. |3TT|) . Thus, both spec- 
troscopy and photometry are consistent with an extincted 
photosphere of a single low-mass star. For object C (oc- 
tagons), however, it is clear that extinction alone cannot 
explain the H — K colour. Object C appears to be 6 mag 
fainter in the K-band image than expected from the pho- 
tospheric fluxes derived from the spectrum. Scattered light 
models provide a m uch better fit. The soli d line shows the 
model MRN12 from iKenvon et"aLl l| 1993^ 1 plus an extinc- 
tion of Av ~ 15 mag, which allows to reproduce the K-band 
flux and the colour simultaneously. This analysis conflrms 
that the object is mostly seen in scattered light through an 
edge-on disk. 

In summary, the information from photometry and 
spectroscopy indicates that AB and C are both low-mass 
stars. Since the mass of AB is likely to be higher than the 
mass of C and AB is located closer to the core of the neb- 
ula, AB is probably the most relevant center of infall in the 
IRAS04325 system. 



~ 50 deg (Fig. [2)| . Note that the position angle of the beam 
is SE to NW, i.e. perpendicular to the orientation of the 
source. The parameters of the gaussian have large uncer- 
tainties, but it is safe to conclude that the position angle 
is consistent with the orientation of the disk, as inferre d 
from the HST images (PA 30 deg, iHartmann et al.lll999l '). 
The same HST images also constrain the radius of the disk 
to 0'.'2 (- 30 AU), which is not resolved by the SMA. The 
elongated structure seen in the submm data might be caused 
by a cold outer disk or by an elongated circumstellar enve- 
lope with a diameter of ~ 1 — 2", i.e. 150-300 AU. This size 
matches well with the constraint from the SED modeling 
(Sect. [521). 

The southern AB source contains a compact component 
that is well matched by a point source, likely caused by a 
small-scale disk. In addition there is a contribution of spa- 
tially extended emission, mostly oriented in E-W direction 
over an area of 3-5" . 



3.4 PSF fitting 

It has been speculated that AB mig ht be a binary (he n ce the 
name). Based on the HST images, IHartmann et al.l (|l999l ) 
argue that 'there is a double structure, possibly indicating 
two sources or possibly indicating a dark absorption lane 
running roughly east- west across the object.' In the HST 
images the separation between the two components is in the 
range of 0'.'2 roughly in north-south direction, corresponding 
to ~ 30 AU. The 'dark lane' would indicate an absorption 
feature in front of the object, possibly a disk with diameter 
of a few tens AU. 

Based on our K-band image, we test for the presence of 
a second source in AB by constructing a model PSF from 
the three well-detected fleld stars outside IRAS04325. All 
sources in the image are fit with this model PSF, using 
daophot within IRAF. The x^ of the fit is < 2 for all field 
stars, 3.2 for object C, and 12.0 for object AB. The contour 
plot of the residuals (Fig. [7]) after subtracting one PSF does 
not show any evidence for a second point source, and fitting 
two PSFs instead of one does not improve the fit. The high 
X^ is most likely caused by the strongly uneven background 
and not by a stellar companion. 

Thus, we prefer to interpret the double structure seen 
in the HST image as an indication for a disk seen at high 
inclination that bisects the image of the star, rather than 
the presence of a resolved companion. This finding is sup- 
ported by the combined information from photometry and 
spectroscopy; for a binary we would expect a mismatch (i.e. 
a later spectral type than expected from photometry). The 
orientation of the 'lane' is roughly perpendicular to the out- 
fiow emanating from AB (Sect.[4|. 

In the submm continuum image both AB and C are 
spatially resolved. The spatial structure is studied after dis- 
carding the data from baselines shorter than 44 m, which 
effectively removes the structures larger than a few arcsec. 
We used the task uvf it in the Miriad package to fit the visi- 
bilities of the sources. Various combinations of point sources 
and gaussian sources were tried for both objects. 

The northern source C is well matched by an elliptical 
gaussian with axes of l'.'2 and 0'.'2 and a position angle of 



3.5 Astrometry 

Based on the PSF fit discussed in Sect. 13.41 the relative 
positions of AB and C were determined with an accuracy 
in the range of O'.'l (1 pixel). The distance between the 
two sources measured from the peak position of the model 
PSF is 8" 32. Simply measuring the positions of the emission 
peaks also gives a consistent distance of 8" 3. From the co- 
ordinates given bv IHartmann et al.l l|l999f ). based on images 
from November 1997, we infer a distance between AB and C 
of 8" 23. Within the errorbar both measurements are consis- 
tent |j The position angles from the available near-infrared 
images agree well: We measure 351.8 deg (E of N), while 
the literat ure values are 351.4 (Connellev et al.ll2008r ) and 
351.6 deg (|Hartmann et al.lll999 ). all with uncertainties of 
;$ 1 deg. 

These measurements constrain the relative proper 
motion between AB and C to < 9masyr~^. For 
comparison, the average proper motion of young stars 
within 5 deg of IRAS04325 is 23masyr~'^ (|Ducourant et al.l 
120051 : iLuhman et al.l 120091 ) , with a standard deviation of 
lOmasyr"^. While the constraints on the proper motion for 
AB and C do not prove that the object is a physically bound 
system, they do confirm common membership in the Taurus 
association. 

The angular distance between AB and C correspond to 
a separation of 1250 AU and a relative movement < 15 AU 
over 11 yr, which results in an upper velocity limit of 6 kms~^ 
between AB and C. This does not rule out the possibility 
that C has been much closer to AB in the past: With a 
velocity of 6kms~^ it would have taken only 10^ yr for C to 
move to its current position from a starting point close to 
AB. It is thus conceivable that C has been ejected by a close 
dynamical encounter with AB in the early evolution of the 
system. 



^ For co mpleteness, the r elativ e distance from the astrometry 
given by IConnellev et al.l (|2008|) based on imaging from 2003-5 
is 8'.'03; the accuracy of their astrometry is not clear. 
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4 MORPHOLOGY OF IRAS04325 

4.1 Structure within IRAS04325 

The IRAS source is located within the small dark cloud 
L1535, a part of the B18 cloud complex in Taurus. The few 
field stars visible in both H- and K-band have H — K colours 
of 2-3 mag, corresponding to extinctions of 25-35 mag, which 
indicates substantial amounts of gas and dust along the line- 
of-sight. 

The near- and mid-infrared images of IRAS04325 
show a bright e longated emission nebula (|Hodapd 1 19941 : 
iHartmann et al.lll99 9: Connelley ct al. 2007]). In our images 
in the JHK bands the dimensions of this nebulosity are 
about 10" X 20", corresponding to 1500 x 3000 AU, with a 
position angle (PA) of ~ 15 — 20 deg. The nebula is divided in 
two lobes. While the northern lobe is brighter, the southern 
one is more extended with two long light streaks visible up 
to 15" to the south of the nebulosity center. The structure 
is best explained as light from object AB scattered by dust 
at the edges of an outflow cavity; indeed, the orientation fits 
with various outflow related features (Sect. [4^ . 

Apart from the two compact sources, the multi- 
wavelength images reveal a complex substructure in 
IRAS04325, see the contour plot in Fig. [T] where we have 
subtracted a point source at the position of AB (see Sect. 
13. 4|) . We note that the PSF subtraction affects only the area 
within 0'.'5 radius of AB and does not alter the structure 
seen in the figure. There are two extended features in close 
proximity to AB (Fig. [3 left panel). The southern feature, 
called D in the following, has a peak 0'.'7 SW of AB and 
appears in the K-band as a half-ring broadly covering the 
south part of the point source. This is best explained by scat- 
tered light from an envelope around the central source, with 
a radius of about 100 AU. The position and structure of this 
feature agrees well with extended CO emission seen in the 
submm regime at radial velocities of 6.5-8 kms~^, indicating 
the presence of cold gas. In the CO map the maximum in- 
tensity is seen at 6.9 kms^^ with a peak at O'.'S south of AB 
(Fig.©. 

The northern feature, called E in the following, appears 
in the K-band as an open cone with an apex located l'.'2 
north to AB, corresponding to 180 AU (see Fig. [T]and[71). 
Two curved light streaks extend this structure towards the 
east and west. The same extensions are seen in the submm 
continuum image (Fig. [21). The position of feature E also 
coincides with a CO lobe seen in a velocity regime between 
2 and 5kms~^ (see Fig. [3}. It is most pronounced at 4- 
4.5kms~^, with a peak position at about I'.'O north of AB. 
Feature E is orientated at a position angle of ~ 15 deg with 
respect to AB, coinciding with the orientation of the emis- 
sion nebula. The morphology and orientation fit well with 
an interpretation as part of the outflow cavity. 

Our own K-ban d image and the one from HST 
IjHartmann et al.ll 19991 ) show a faint S-shaped structure con- 
necting the point sources AB and C and extending further to 
the north (Fig. [T]and[7l right panel). This is most likely scat- 
tered light from dust features. It is possible that this 'bridge' 
between AB and C is coincidental; we see the extension of 
the outflow-related feature E (see above) and light streaks 
originati ng from source C in su perposition in the plane of 
the sky. IHartmann et al.l l| 19991 ) argued that the two faint 
features near C are scattered light from the cavity edges of 





Figure 8. Spitzer IRAC images centered on the nebulosity north 
of IRAS04325 at 4.5 (left) and 8.0 /im (right). IRAS04325 is 
marked with a black circle. Noteworthy are a) the emission knot 
seen at 4.5 /im (left panel, marked with arrow) at the northern 
end of the nebula and b) the dark feature visible at 8 ^m. 



an outflow from object C. Their direction is 60-90 deg offset 
from the orientation of the disk seen in the HST images, 
which flts into the outflow interpretation. However, in our 
CO observations there is no evidence for outflow lobes at 
the same positions (see below). 

Therefore we prefer a scenario in which these structures 
belong to a physical 'bridge' between AB and C. Such a 
'bridge' could be indicative of gravitational interaction in 
the early evolution of the system, but the long-term stabil- 
ity of the outflows (Sect. 14. 2p argues against this interpre- 
tation. Instead we propose to explain the structure as dust 
fea tures from interacting accretion streams (as speculated 
bv IHartmann et al.l (jl999f )). The S-shaped morphology in- 
deed resembles the dens ity maps in accr e ting b inaries from 
the SPH simulations bv lBate fc Bonnelll (|l997l '). 

Two CO lobes are detected in close proximity to source 
C in a narrow velocity range of 6.5-7.5 kms~^ (see Fig. |3]). 
The strongest one is almost exactly coinciding with the posi- 
tion of the compact source seen in the continuum image: the 
offset between peak of the CO intensity and peak in the con- 
tinuum is < I'.'O. The PA for the elongation of this CO lobe 
is roughly 300 deg. The second CO blob is located roughly 
5" away at PA ~ 290 deg with respect to C. A third CO de- 
tection is seen in the velocity range of 4.4-4.8 kms~^ 5" west 
of object C at PA of ~ 270 deg. The orientation of these 
three CO features is roughly perpendicular to the alignment 
of the disk (PA 30 deg). Thus, they could be associated with 
an outflow from source C. 



4.2 Structure outside IRAS04325 

A region of diffuse near-infrared emission is found about 
40" north and 15" east of IRAS04325AB, corresponding to 
a PA of 20 deg (see Fig. |8} . Thus, this area is located in 
the extension of the axis defined by the elongated emission 
nebula. The feature h as been mentioned befo re in the liter- 
ature, for example bv IPark fc KenvonI (|2002l ). and is most 
pronounced at 2 to 5/im. In our K-band image the feature 
has a size of 20" x 20" . The proximity to IRAS04325 and the 
location along the outfiow axis lead us to believe that the 
feature is physically associated with IRAS04325. Under this 
assumption the physical distance from IRAS04325 would be 
~ 6000 AU. This patch of emission could be scattered light 
at the northern end of the outfiow cavity. 
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Figure 7. Contour maps of the K-band image from Gemini/NIRI around tlie source IRAS04325. One small tickmark corresponds to 0'.'6 
(5 pixel). The left panel shows a small area around source AB after subtracting one point source centered at the peak of the emission. 
The total image size is 6" X 6" . The two extended structures D and E to the south and north of the point source are clearly visible. 
The right panel shows the northern part of the source around object C; image size 10" X 10". The dashed lines mark the light streaks 
connecting AB and C and extending further to the north. Note that the dark lane from the disk around C, as seen in the HST image, is 
not resolved here. 



In the IRAC bands at 3.6 and 4.5 ^m a 10-15" wide, 
dark band is visible between the IRAS04325 nebula and the 
emission region further north (Fig. |8l left panel), possibly a 
foreground cloud blocking the view to the source. 

At 8/^m the bright diffuse patch in the north of 
IRAS04325 is not visible, instead the image shows a roughly 
circular, 1.5' x 1.5' absorption feature centered 70" north 
and 30" east of IRAS04325 (PA 23 deg), i.e. just north of 
the bright patch (Fig.[8l right panel). This dark feature has 
fainter extensions to the east and south, overlapping with 
IRAS04325. Its center is located along the outflow axis of 
IRAS04325AB. The position agrees weU with the core L1535 
N- SMM seen in the submm map s of IRAS04325 presented 
bv iHogerheiide fc Sandelll (|2000|), indicating that this fea- 
ture represents an overdensity of dust. It is possible that 
the northern lobe of the outflow is blocked by this core, ex- 
plaining the lack of large-scale outflow features further to 
the NE (see below). 

The dark feature is only visible at 8 fim, and not at any 
other wavelength in the mid-infrared. In the 8 fim image it is 
the only such feature within a radius of one degree. This may 
be related to the presence of strong PAH e mission bands in 
the other IRAC bands (e.g. IPeeters et al.1.2002 .'). It is con- 
ceivable that we see a region of strong absorption projected 
against the bright PAH back groun d, as they are often found 
in the Galactic plane (e.g. .E gan et al. 1998). In fact, the 
combination of submm emission and 8fim absorption has 
been observed for a nu mber of cores in infrared dark clouds 
jVasvunina et al.ll2009h . 

IRAS04325-I-24 02 has long been known to drive 
molecular outflows (iHeyer et al.1 Il987l : iTerebev et al.l 1 19891 : 
iMoriartv-Schieven et al.lll992l V Speciflcally, a large body of 
evidence is available for an outflow driven by source AB at 
position angle of 10-20 deg. Features associated with this 
outflow are the elongated emission nebula (also seen in 



HCO^. IHogerheiide et al.lll998l ). the scattered hght feature 
E (Sect.[4T]), the HH objects 434-436 about 1 deg south- 
west (2.7 pc) of the IRAS source at PA of 17.5-19.5 deg 
jWanget al.ll200ll ). and HH703 30' so uth of IRAS04325 
(1.4 pc) at PA of 16 deg (|Sun et al.ll20o3 ') 

We searched in the H2 images as well as in the IRAC 
images for more shock features along the outflow axis. In 
particular the IRAC band 2 (4.5 fim) is known to con- 
tain a number of bri ght emission lines of shocked H2 (see 
lYbarra fc Ladall2009l . and references therein). The group of 
Herbig-Haro objects 434-436, so far only detected in the [SI I] 
line and possibly part of a bow shock (|Wang et al.ll200ll ). 
is well-visible in H2 and in all four IRAC bands. In the 
4.5 fim image we flnd an additional knot about 1.1' north of 
AB (10000 AU) at PA 14 deg (a 4''35'"36':5, 6 -K24°09'25'.'5 
(J2000), Fig. IS] left panel). This feature is clearly extended 
and not seen at any other wavelength, which argues for an 
association with the outflow. No other obvious feature was 
found for PA of 15-20 deg and separation from IRAS04325 
up to 1 deg in both directions. 

There is a second large-scal e outflow em a natin g from 
IRAS04325-h2402, identifled by iHever et al.1 (|l987D . The 
outflow is seen in integrated CO emission spanning 12', 
i.e. 0.5 pc if related to IRAS04325. It is described as well- 
coUimated, redshifted, and monopolar. In clear disagreement 
to the first outfiow discussed above, this second flow is ori- 
ented from the IRAS source towards the NW at position 
angle around ~ 310 deg, roughly perpendicular to the disk 
of source C. 

The most likely source for this second outflow is object 
C. This star clearly has jet activity, proven by the variety 
of outflow related emission lines in the near-infrared spec- 
trum. Three CO lobes are seen around this object with PA 
of 270-300 deg, rough ly matching the orientation reported 
bv iHever et al.l ([1987|), which could belong to the outflow 
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as well. We do not find any additional counterparts of this 
outflow in the H2 and the IRAC images. 



5 SED MODELING 

We aim to reproduce the spectral energy distributions 
(SED) for objects AB and C with a Monte Carlo radiative 
transfer code . In ad dition to the data presented here and in 
IScholz et al.l (12003) we make use of the Spitzer /IRS spec- 
trum available for object AB (JFurlan et al.ll200a ). Only the 
SL exposure is used, which is obtained through a 3'.'6 wide 
slit oriented roughly in N-S direction and covers the wave- 
length range from 5.2 to 14 /^m. This spectrum matches well 
with the photomery at 11.2 /iin. We do not use the LL spec- 
trum (obtained through a 10'.' 6 slit) and the IRAC/MIPS 
photometry, due to their poor spatial resolution. 

Our model code computes thermal radiation and non- 
spherical scattering from dust in a two-dimensional axisym- 
metric system with central source, flared disk, and larger- 
scale flattened envel ope with a narrow ev acuated bipolar 
outflow cavity (e.g. IWhitnev et alj l2003bl lal). The central 
photospher es are approximated with NextGen mod el at- 
mospheres iHauschildt et al.lll999l : lAllard et al.ll200ll ') with 
log (7 = 4. The main free parameters for the disk are mass, 
outer radius, inner radius, accretion rate, scale height fac- 
tor ho and flaring power (5. The latter two are defined as 
h — ho{r/Ri,)^ , with h the scaleheight of the disk and r the 
distance from the star. The parameters for the envelope in- 
clude envelope accretion rate, mass, outer radius, and out- 
flow cavity angle. In addition the model fits the ambient 
density and the interstellar extinction between the source 
and the observer. 

For the envelope we assume dust w ith an interstellar - 
like size distribution, using the model bv iKim et al.l (|l994h . 
In the disk we use a dust model including silicates and car- 
bonaceous grains with solar abundances. Here the grain size 
distribution is a power law with an exponential decay for 
particles with sizes above 50 /xm and a formal maximum 
grain size of 1 mm (|Wood et al.l l2002h . This dust model 
thus allows for the presence of larger grains than in the 
ISM. For more d etails and ref e rence s abo ut the radiation 
transf er code see IScholz et al.l (|200y) and iRobitaille et al.l 
(|200g). The code has recently been extended to treat the 
non-thermal reprocessin g of photons by PA H molecules and 
very small grains fVSG. IWood et al.luOOSi '). usi ng opacities, 
abund ances, and re-emission templates from iDraine fc Lil 
(|2007l ). 

Given the large number of free parameters and the nu- 
merous simplifying assumptions (e.g., two-dimensional disk, 
axisymmetric system, n o variability), the modeling is sub- 
ject to degeneracies fe.g. lChiang et al.ll200ll : IRobitaille et al.l 
[2007). The models presented here are by-eye fits to the SED 
and do not represent a unique solution. We do not claim to 
have explored the full parameter space. Instead, we use the 
model to verify if the SEDs can in principle be understood 
in the commonly used framework of disk/envelope models. 
In addition, this will provide an independent test for the 
plausability of the constraints for the physical parameters 
of the system. The parameter best constrained by the mod- 
eling are the circumstellar mass, as it depends only on the 
submm/mm datapoints. On the other hand, the mass dis- 
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Figure 9. SED fitting for IRAS04325AB. Shown with sohd line 
is a model calculated for a square aperture of 200 AU, i.e. compa- 
rable to the observations. The dash-dotted line shows the same 
model for the full scale of the simulation (10^ AU). The photo- 
spheric spectrum is plotted as dashed line. The observed fluxes 
are plotted with triangles, except for the IRS spectrum, which 
is seen as a thick solid line between 5 and 14 ^m. The models 
assume an external (i.e. outside disk and envelope) extinction of 
Ay = 3 mag. 



tribution in disk and envelope, the disk inclination, as well 
as the extinction are poorly constrained. 



5.1 Object AB 

Figure |9] shows a model for object AB computed with- 
out the inclusion of PAH/VSGs. We find that the two- 
dimensional model with disk, envelope, and outflow cavity 
is able to reproduce many of the features of the observed 
SED. One matching SED model is shown in Fig. [21 This spe- 
cific model is based on the following parameters, consistent 
with the observational constraints: stellar mass 1 Mq , stellar 
radius 2R0, effective temperature 3800 K, disk inclination 
78deg, disk mass lO^'^M©, disk radius 50 AU, disk accre- 



tion rate 10 * Mr yr ^, scale height factor 0.01 R*, flaring 



power 1.25, envelope accretion rate 3 ■ 10 ^ M© yr ^ , cavity 
half opening angle lOdeg. The total circumstellar mass is 
0.18 Mq for the fuU scale of the simulation of lO" AU. 

These results are consistent with previous a ttemp ts to 
model the source. For example, IWhitnev et al.l (| 19971 ) find 
a centrifugal disk radius of 50 AU, an inclination of 72- 
90deg, and a total infall rate of 5 • 10~^ M© yr~^ based on 
near-infrared imagi ng polarimetry. The model presented by 
iFurlan et al.l l|2008l l aims to reproduce the full IRS spec- 
trum and broadband photometry from the near-infrared to 
the mm regime and results in a centrifugal disk radius of 
100 AU, a cavity semiopening angle of 15 deg, and a disk 
inclination angle of 80 deg. 

The largest discrepancies between model and observed 
SED occur in the near- and mid-IR. This part of the SED 
is sensitive to small changes in extinction, stellar parame- 
ters, disk inclination an d to non-axisymmetric geometries 
(jlndebetouw et aLll2006l ). The images of this object show a 
complex geometry, i.e. 3D effects will likely affect the mid-IR 
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Figure 10. Modeling tiie IRS spectrum (solid line) for 
IRAS04325AB. Tlie dotted line is a model without PAH, i.e. the 
same as in Fig. |9] For the dash-dotted line model, wo include 
PAH in the disk. The dashed line has PAH in the disk and en- 
velope. Triangles show the photometric datapoints. All models 
are calculated using the same parameters as in Fig. |9] They as- 
sume an external (i.e. outside disk and envelope) extinction of 
Ay = 3 mag. 



The IRS spectrum for this object, shown in detail in 
Fig. 1101 is not well-reproduced by our model. The deep ab- 
sorption trough around 9/im is most likely caused by the 
amorphous silicate abs orption band at 9 .7 /im, as seen in 
many other protostars (JFurlan et al.ll2008l l. In addition, the 
spectrum could be affected by the frequ ently observed ice ab - 
sorption bands at 6.0 and 6.8 /im (e.g. Ivan Dishoeckl 12004 ). 

On the other hand, the peaks at 6.2 and 7.6 fim coincide 
with tw o of the most prom inent features of PAH emission 
spectra (jPeeters et al.ll2002l ). We made an attempt to model 
the spectrum by including treatment of PAH/VSGs in the 
radiative transfer simulations. In Fig.[TO]we show the model 
with PAH only in the disk (dash-dotted line) as well as PAH 
in disk and envelope (dashed line). Qualitatively the shape 
of the two peaks is matched by the model. The absolute 
flux levels, however, are not adequately reproduced by these 
PAH models. The strength of the peaks is highly sensitive 
to the location of the PAH; in this case, a model with PAH 
only in the disk provides the best match. 

PAH emission is generally found to be rare in disks 
aroun d young stellar objects, particularly in late-type ob- 
jects (|Geers et al.ll2006l : Furlan et al.1 120061 ) . with detection 
rates below 10%. Given the rarity of PAH features and the 
ubiquity of ice features in low-mass protostars, we consider 
it more likely that ices are predominantly responsible for 
the shape of the mid-infrared spectrum in IRAS04325AB. 
As our analysis has shown, the presence of PAH provides an 
alternative plausible interpretation for the peaks. To disen- 
tangle these effects, it is nec essary to include ice opacities in 
the models, see for example iGrapsi et al.l (|2008r ). 



5.2 Object C 

In Fig. [TT]we show the SED for object C and a model SED 
that matches most of the datapoints. In the modeling we fix 
the effective temperature of the central object (3400 K, see 



Figure 11. SED fitting for IRAS04325C. Shown is a model 
with Ay = 17.5 mag that matches all available independent con- 
straints (effective temperature, disk radius, accretion rate). Note 
that the J-band flux is an upper limit. The dashed line shows the 
photospheric spectrum, the dotted line the intrinsic SED without 
external extinction. 



Sect.|3]2l), the accretion rate (10~'* Mayr~\ see Sect.lT2l) . 
and the outer disk radius (30 AU, iHartmann et al.lll999l ). 
In addition, we require a high disk inclination (i > 80deg). 
An edge-on or close to edge-on geometry is necessary to 
reproduce the appearance of the HST image and to explain 
why photometry and spectrosocpy cannot be simultaneously 
reproduced by an extincted photosphere (see Sect. 13.3")) . 

The model shown in Fig. [11] is calculated for i?* = 
1.8 R0, i = 85deg, ho = O.OIR*, /? = 1.25, i.e. the disk 
is assumed to be affected by relatively strong flaring. The 
total mass in disk and envelope has to be between 0.001 and 
0.01 M0 to reproduce the submm/mm datapoints (assuming 
the canonical gas-to-dust ratio of 100). The model predicts 
most of the circumstellar mass to be in the disk (0.004 Mq), 
with only a thin envelope (6 ■ 10~^ M0), but other conflgu- 
rations are conceivable. The envelope has an outer radius of 
300 AU, an accretion rate of 10~^ M© yr~^, and an outflow 
cavity with half-opening angle of 15 deg. 

The model assumes interstellar extinction of Av = 
17.5 mag. This does not include the extinction caused by 
disk and envelope, which is extremely high for this object 
due to the edge-on disk. The value for the extinction agrees 
well with the estimate we obtain from the combined infor- 
mation of near-infrared spectroscopy and photometry (Sect. 
I3.3|l . It differs signiflcantly from the extinction used in the 
SED modeling for object AB {Ay = 3 mag. Sect. 15.1) . which 
is implausible because outside the disk/envelope systems the 
extinction should be comparable. As noted before, the con- 
straint for Av obtained from SED modeling is subject to 
large uncertainties, due to the numerous degeneracies in the 
near-infrared SED. 

As seen in Fig. 1111 the photometric datapoint in the 
M'-band at 4.7 /im is difficult to reproduce with these 
models. This does not depend on the particular choice of 
disk/envelope properties; the datapoint is an outlier for all 
reasonable combinations of parameters. Most models predict 
a decreasing SED (in XFx) between 3 and 5/im, whereas the 
datapoint at 4.7 /im is 30% higher than the one at 3.6 /im. 
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One possible explanation is additional flux generated by the 
outflow, which is not part of the model. The M'-band con- 
tains a number of strong H2 emission l ines typical for shocks 
in protostellar outflows (|Ybarra fc Lada 200 9). in particular 
the 0-0 S(9) transition at 4.69 /im. Shocked H2 emission, as 
seen in the near- infrared for this object (Sect. 13. 2p . could 
contribute significantly to the M'-band flux. 



6 SUMMARY AND CONCLUSIONS 

We present a multi-wavelength study for the embedded pro- 
tostellar system IRAS04325+2403. IRAS04325 is one of very 
few known multiple YSOs with at least one resolved disk and 
two protostellar jets. New imaging with subarcsecond reso- 
lution has been obtained in the infrared covering the wave- 
length range 1-12 /^m and in the sub- millimetre at 870 /im. 
This is combined with low-resolution near-infrared spec- 
troscopy for the two compact sources in the system and 
with literature data. In the following we summarize the con- 
straints on the properties of IRAS04325: 

(i) We show that the both components in the system are 
low-mass stars. Object C is an M dwarf, probably with an 
effective temperature around 3400 K, while AB is likely of 
earlier type and more massive. Object C is unlikely to be a 
brown dwarf, as speculated previously. 

(ii) Despite the wealth of available data, there are enor- 
mous difficulties in estimating the fundamental parameters 
of the central sources. In particular, the presence of veiling, 
disks, and strong line-of-sight extinction puts limits on our 
ability to reliably derive effective temperatures from spec- 
tra. The poorly constrained age hampers the conversion to 
masses. 

(iii) The near-infrared emission line spectrum for object 
C provides strong evidence for ongoing accretion and jet 
activity. 

(iv) We argue that the previously seen double structure 
of object AB is not due to the presence of a close stellar 
companion. Instead it is best explained by an absorption 
lane in front of the object, most likely caused by a disk. 

(v) Com paring our astrometry with the positions re- 
ported by iHartmann et al.l (|l999l ). we put a limit of < 
9masyr~^ on the relative proper motion between AB and 
C. This confirms common membership in the Taurus star 
forming region. 

(vi) The morphology and SED for object AB is consistent 
with the presence of a small disk (radii ~ 50 AU radius) and 
a larger envelope (< 100 AU radius). 

(vii) Object C has been found to have an edge-on disk 
with radius 30 AU by IHartmann et al.l (|l999f ) . Our obser- 
vations confirm that the object is mostly seen in scattered 
light due to the high inclination. Based on submm imaging 
and SED analysis we identify an elongated envelope or outer 
disk with radius of < 300 AU around this object. 

(viii) The mass of circumstellar material in the system is 
well-constrained by the submm/mm datapoints. While the 
disk masses are quite small (0.001-0.01 Mq), the total mass 
within 10000 AU is in the range of 0.2 M©. 

(ix) Many of the features in the SED and the multi- 
wavelenght images can be understood in the framework of 
commonly used two-dimensional disk and envelope models. 
The remaining inconsistencies may well be caused by the 



3D structure of the system, which is visible in the images 
but not accounted for in the models. The mid-infrared spec- 
trum for source AB shows spectral features that cannot be 
matched with a standard model using a single population of 
dust grains; instead they require the inclusion of silicates, 
ices, and possibly PAHs. 

(x) We identify several new features associated with a 
parsec-scale outfiow driven by object AB with position angle 
of 10-20 deg. Object C likely has its own large-scale outfiow 
at position angle of 270-310 deg. 

(xi) Various pieces of evidence indicate an early evolu- 
tionary state for IRAS04325. The star-disk systems are em- 
bedded in a core with a large line-of-sight extinction of 
Av ~ 30 mag. The object is und etected in the Taurus XMM 
X-ray survey (jGudel et al.ll2007l ). which is unusual for Class 
II objects (15% undetected), but common for protostars 
(60%). The bolometric luminosity of t he system, as de ter- 
mined by fitting the unresolved SED (|Froebrichll2005h . is 
ar ound 1.0 Lp. According to the pre-main sequence tracks 
bv iD Antona fc Mazzitellil (|l994l . update 1998), two objects 
with effective temperatures around 3400 K would have a 
total luminosity of 1.3, 1.0, 0.91.0 at 10^, 2 x 10^, and 
5 X 10"" yr. Thus, an age of a few times 10^ yr seems rea- 
sonable. 

(xii) The position angles of disks, envelopes, and outflows 
put useful limits on the geometry of the system. For the AB 
disk/envelope, C disk/envelope, and binary orbit the po- 
sition angles are approximately 90, 30, and 350 deg. Thus, 
the disk/envelope systems in IRAS04325 are strongly mis- 
aligned with respect to each other (60 deg) as well as with 
respect to the orbital plane of the binary (80 and 40 deg). 
The disks are neither aligned nor coplanar. 



The last flnding, flrst stated by IHartmann et al.l (|l999r ) 
and bolstered here by more evidence, warrants some further 
comments. From polarimetry studies it is known that the 
dominant majority (80-90%) of young binaries has disk ori- 
entat i ons within 30 deg of the orbital plane (|jensen et al.l 
|2004 iMonin et all I2OO6I). There a re exceptions, found by 
imag ing the le ts (e.g. iDavis et al. 1994; Lim & Takakuwa] 
I2OO61 ) or disks (iKang et al.ll2008l ). but many of them seem 



to be higher-order systems where one o f the components is 
again a close binary (see the review bv lMonin et al.ll2007l l. 
This could be the case for IRAS04325 as well, although we 
argue against the presence of an additional component. 

Misalignments can be either primordial or arise dur- 
ing the early evolution, due to dynamical interactions in 
a multiple system. In the case of IRAS04325 the presence 
of parsec-scale outflows with constant directions to within a 
few degrees seems to exclude strong alignment changes in the 
recent history. Typ ical dynamical timescales for parsec-scale 
jets are lO'^yr Ce.g. lEisloffel fc Mundtl^l997^ . For the NS jet 
driven by object AB we estimate a lifetime of lO"* — 10^ yr, 
assuming typical velocities on the order of a few lOOkms"^. 
The dynamical timesc ale for the NW je t has been deter- 
mined to be 2 X lO'^yr (|Hever et al.lll987h . Thus the system 
has been stable over a substantial fraction of its early history. 
This is confirmed by the absence of strong relative motion 
between the two components, as it might be expected after 
recent gravitational encounters. 

In general the misalignments between disks and or- 
bits in binary systems will decay with time; the timescale 
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for this readjustment is expected to be i n the order of 20 
orbital periods or 10^ yr for IRAS04325 (JBate et all I2OO0I : 
I Jensen et al.ll2004l l. which might be longer than its lifetime. 
For higher-order systems the timescales for realignment are 
much longer, because the various orbital planes have to be 
realigned, not only the disks. In addition, the disk radii in 
IRAS04325 are significantly smaller than the orbital separa- 
tion (< 50 AU vs. 1250 AU). In such cases it is possible that 
the tilt between disks and orbital plane grows with time 
and an initial misalignment is enhanced l|Lubow fc Qgilvid 
bOOd ). Thus, IRAS04325 might be a good case for a system 
where we see primordial misalignment between disks and 
orbital plane. 

For a wide, isolated system like IRAS04325-)-2402 there 
are two main possible avenues of formation: fragmentation 
from a centrifugally supported core (rotational fragmenta- 
tio n) or from turbu l ent m aterial (turbulent fragmentation) - 
see lGoodwin et al. (2007) for a review. As discussed exten- 
sively bv iBate et al., (,2000. ). primordial misalignments are 
easy to produce in a turbulent fragmentation scenario, for 
example due to spatial variations in the direction of angular 
momentum within the core. On the other hand, rotational 
fragmentation predicts a strong alignment between disks 
and orbit. The picture presented here for IRAS04325 shows 
that strong misalignments can be present even in small-scale 
star forming environments and favours turbulent fragmen- 
tation as the formation scenario. 
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